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Abstract 

An Extended Kalman Filter is developed to estimate the linearized direct and indirect stiffness and 
damping force coefficients for bearings in rotordynamic applications from noisy measurements of the 
shaft displacement in response to imbalance and impact excitation. The bearing properties are modeled as 
stochastic random variables using a Gauss-Markov model. Noise terms are introduced into the system 
model to account for all of the estimation error, including modeling errors and uncertainties and the 
propagation of measurement errors into the parameter estimates. The system model contains two user- 
defined parameters that can be tuned to improve the filter’s performance; these parameters correspond to 
the covariance of the system and measurement noise variables. The filter is also strongly influenced by 
the initial values of the states and the error covariance matrix. The filter is demonstrated using 
numerically simulated data for a rotor-bearing system with two identical bearings, which reduces the 
number of u nk nown linear dynamic coefficients to eight. The filter estimates for the direct damping 
coefficients and all four stiffness coefficients correlated well with actual values, whereas the estimates for 
the cross-coupled damping coefficients were the least accurate. 

Introduction 

Compliant surface foil air bearings are self-acting hydrodynamic bearings that support a rotating shaft 
on a thin layer of gas between the shaft surface and a sheet metal foil supported by a series of stiff bump 
foils (see fig. 1). Advancements in coating materials and methods to design the supporting bump foil 
system have increased the load bearing capacity in gas foil bearings (DellaCorte and Valeo, 2000) and 
improved their stiffness and damping properties (Salehi et al., 2004) and their ability to operate in high 
temperature environments (DellaCorte et al., 1999; Heshmat et al., 2007). These advancements have 
made gas foil journal bearings prime candidates for implementation in new applications critical to 
NASA’s mission, such as systems to convert nuclear energy to electricity for power generation in space 
and more efficient, lower maintenance propulsion systems for aeronautics. Additionally, gas foil bearings 
are increasingly being considered and used in commercial applications such as microturbine generators 
and industrial blowers. Coupled with these materials and design advancements is the need for 
experimentally verified tools for predicting steady state performance and rotordynamic properties of gas 
foil bearings. The capability to measure the linearized rotordynamic properties of actual gas foil bearings 
is necessary to validate and calibrate theoretically predictive tools and to boost the confidence that high- 
speed rotordynamic machinery utilizing these components will operate effectively and reliably. 
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Figure 1 . — Cross section of a gas foil journal bearing. 

Knowledge of the bearing dynamic properties is instrumental in predicting the performance 
characteristics of high speed turbo machinery. As a shortcut to completely modeling the physics of the 
bearing, the bearings are simply represented by linearized springs and dampers, in which case the 
challenge then becomes to estimate the linearized force coefficients for these components. Bearing 
reaction forces in this form are most convenient for analyses involving the whole rotor-bearing system for 
which each bearing is just one of many integrated components. 

Several discrete-time domain and frequency domain techniques have been presented in the literature 
for estimating the linearized stiffness and damping (and inertia) coefficients for bearings and seals. Tiwari 
et al. (2004) provides a review of the literature on the topic. Discrete-time domain methods, which 
generally utilize least-squares estimators, often suffer from the disadvantages that the bearing estimates 
are biased and highly sensitive to measurement noise and errors introduced by transforming the 
continuous time model into a discrete-time model (Sahinkaya and Burrows, 1984). These problems are 
mostly overcome, resulting in improved bearing parameter predictions, by employing estimation methods 
in the frequency domain. Most of the methods discussed in the literature estimate the u nk nown bearing 
parameters using motion measurements from systems that are excited by some form of known external 
forcing, e.g., shakers or impacts. Experimental systems with the capability for external excitation from 
shakers are costly and generally limited to the laboratory, whereas many rotating machinery in use are 
instrumented with shaft displacement sensors and housing accelerometers for condition monitoring 
purposes. For systems without the capability of external harmonic excitation, estimation methods have 
been developed that rely on imbalance response measurements or the response from impact excitation. 
The system of equations generated by imbalance excitation tends to be ill-conditioned, especially if the 
orbit is near circular (Roberts et al., 1990; Murphy and Wagner, 1991), and the results from these 
methods often show considerable scatter (Tiwari et al., 2004, De Santiago and San Andres, 2007b). 
Bearing identification from impact excitation has demonstrated promising results (Qiu and Tieu, 1997; 
Tiwari et al., 2004, De Santiago and San Andres, 2007a). Also, Howard et al. (2001) determined direct 
bearing dynamic stiffness and damping coefficients for a gas foil journal bearing at the system natural 
frequency by measuring the free vertical displacement of the shaft in response to an impact excitation. 

The system natural frequency provided the stiffness estimate, and the damping was estimated from the 
decay rate. 

The objective of this paper is to demonstrate the application of an Extended Kalman Filter (EKF) to 
estimate eight linearized rotordynamic (stiffness and damping) force coefficients for bearings from 
measurements of the system response to imbalance and impact excitation. The EKF is a time domain 
filter that provides an efficient recursive algorithm for estimating states in a nonlinear dynamical system 
from noisy measurements. In this application, the dynamical system states to be estimated must include 
the bearing stiffness and damping properties, which results in a system governing equation that is a 
nonlinear function of the states. As a consequence, the estimation algorithm will employ the Extended 
Kalman Filter instead of the traditional Kalman Filter, which is used for linear systems. Fritzen and 
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Seibold (1990) have successfully applied the EKF to identify the damping factor and depth parameter for 
a simple rotor with a cracked shaft and the inertia, damping and stiffness coefficients for annular seals in 
a high pressure turbopump test rig. 


Dynamic Model of the Shaft and Bearings 

A schematic of the rotor-bearing system is shown in figure 2. The shaft layout is chosen to match the 
rotordynamic simulator facility at NASA’s Glenn Research Center shown in figure 3 (Howard, 2007). In 
future work, the principles developed in the current paper will be applied and compared to experimental 
data obtained from the rotordynamic simulator facility. The ratio of the shaft radius to length is relatively 
large, so the rotor is assumed to be rigid. The rotor displacement is represented by four coordinates: X\ 
and iq correspond, respectively, to horizontal and vertical displacements from equilibrium (x U) andy 10 , 
respectively) of the shaft at the center of the bearing located near the turbine disk, and X 2 and v 2 are 
similar shaft displacements from equilibrium (x 2 o and >’20, respectively) at the bearing near the thrust 
bearing plate. The linearized bearing model includes only the elastic and dissipative effects; inertial 
effects are ignored. The x and y direction bearing forces (f, x and fb y , respectively) on the shaft are 
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The bearings are assumed to be identical, so they are modeled with the same rotordynamic coefficients, 
which minimizes the number of parameters to be estimated to a total of eight, four direct (k xx , k vv , c xx , and 
c yy ) and four cross-coupled (k xv , k vx , c xy , and c yx ) coefficients. 

The linearized equations of motion for the shaft can be expressed as 

Mz + (C + QG)z + Kz=F 0 (t), (2) 


where z = {xi,yi,x 2 ,y 2 } T » z(t = 0) = z 0 , and z(t = 0) = z 0 . The mass and gyroscopic matrices, M and 
G, respectively, are 
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where L = l x + l 2 . Likewise, the stiffness and damping matrices, S and C, respectively, are 
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Figure 3. — Schematic of the NASA Glenn Research Center rotordynamic simulator test rig. 


The forcing terms in the equations of motion result from the effects of imbalance and impact 
excitation. Imbalances ni\i\ and m 2 r 2 are oriented at angles of <|)i and (j) 2 , respectively, in planes located at 
d\ and cA, respectively, from the center of mass. Likewise, an external impact force with x and y 
components, respectively, of f x and /,'., acts at a distance of z p from the center of mass toward the thrust 
bearing plate, leading to 
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The values for the imbalance represent both the residual imbalance of the shaft and any applied 
imbalance. The residual imbalance can be measured with a balancing machine, and the applied imbalance 
is assumed to be known. The impact force is also assumed to be known as a function of time, which can 
be achieved if the impact is performed with an instrumented impact hammer. 

In most situations it is not possible to place the displacement transducers in a location so that they 
measure the shaft motion at the centers of the bearings. In the case when the displacement transducers are 
located at p\ and p 2 from the shaft center of mass, see figure 2, the shaft motion at the bearing centers, 
denoted by x\,y\, x 2 , and y 2 , can be computed from the transducer measurements, x p \,y p \, x p2 , and y p2 , with 
the following linear transformation, 
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For the EKF, the bearing rotordynamic coefficients must be included as states so they can be 
estimated; consequently, the global state vector becomes 
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where r contains the rotordynamic coefficients, such that r = j c xx ,c xy ,c yx ,c yy ,k xx ,k xy ,k yx ,k yv ^ . Within 

the framework of the EKF estimation scheme, all of the variables in the state vector, including the 
rotordynamic properties, will be modeled stochastically. The rotordynamic coefficients are assumed to be 
constant, but, realistically, the estimates will vary with time due to some random error within the 
measurements. The error within the estimates will be parameterized by the array of stochastic variables, 
w r . A first-order Gauss-Markov model for the rotordynamic parameters is 

f = w r . (8) 

The vector w r is assumed to be Gaussian white noise with zero mean and covariance matrix Q r . 

Now, the equations of motion are combined with the Gauss-Markov models to form a compact, 
nonlinear system model, 
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Here, w v and w z also represent process noise vectors, assumed to be Gaussian white noise, with zero 
means and covariance matrices Q v and Q z , respectively. The process noise vectors form 

with a corresponding covariance matrix Q. These noise terms account for all of the 


T 

W = 


T T T 

w v w z w r 


error within the estimates, including modeling errors and uncertainties and the propagation of 
measurement errors into the estimated values. The system model can be represented as a discrete-time 
equation for discrete time steps t k (k = 0, 1,2 . . .) by 


S A+ 1 ~ s k 


+ f fi+I /(s,F 0 ,f)<# + w Jfc . 

J h 


( 10 ) 


Finally, an output vector, y k , is formed from discrete samples of the shaft displacement. These 
measurements will also be contaminated with noise, so the model for the output vector is comprised of a 
linear combination of the true shaft displacement states, s z / f , and a stochastic random noise parameter, \ k , 
to compensate for the measurement noise, such that 


y,t= s Z £ +v *= Hs /fc + v /t, (ii) 

where H = [0 4x4 I 4x4 0 4x8 ] . The measurement noise parameter, v*, is assumed to be white noise with a 
constant covariance matrix R. 


The Extended Kalman Filter (EKF) 

The EKF is a predictor-corrector type of estimation scheme that provides estimates at each time step 
for the states, , and the error covariance matrix, P /; . The error covariance matrix is the expected value 
of the error between the true states and the state estimates. 

The first step in the process to project the state estimate from one time step to the next is to predict the 
projection based solely on the model and the old state estimate, forming intermediate predictions of the 
state and covariance matrix at the new time step. To distinguish between predictions and estimates at a 
time step, predictions of parameters are designated with the over-bar symbol and estimates are marked by 
an over-caret symbol. The state prediction at the new time step is computed by solving the original 
nonlinear governing equation, giving 


»A+l =h + r* +I /(s,F 0 ,f)<#. (12) 

J h 

Using the “Continuous-Discrete EKF” formulation for nonlinear systems with discrete measurements, as 
described by Gelb (1974), the error covariance prediction is determined by 

P* + i = ^ + j^ +1 (AP + PA r + Q)dt . (13) 

The time integrations in eqs. (12) and (13) are computed using an efficient high-order numerical 
integration scheme. The A matrix in the error covariance prediction represents a linearization of / about 
the current estimate of the state vector, 

\ = (^A’PqaTa) ( 14 ^ 

5s 

Initial values for the state vector and error covariance matrix are set to s 0 and P 0 , respectively. 
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The second step is to update the state and error covariance predictions using the measurement data, 


h+\ - s /t+i +K /t+i(yA+i ~ Hs /t+i) O 5 ) 

and 

P* +1 =(l-K* + iH)P* +1 . (16) 

The term K^+i is known as the Kalman gain factor. The following formulation for the Kalman gain factor 
blends the state estimates with the measured data in a manner that minimizes the error covariance 
(Friedland, 1986), 


K/fc+i = P, + iH r (HP, +1 H r + r)" 1 . (17) 

Equations (12) to (17) are repeated at consecutive time steps corresponding to each discrete measurement, 
providing sequential, filtered estimates of the shaft displacement variables, their derivatives, and the 
rotordynamic coefficients. 

Though the EKF is relatively simple and computationally similar to the traditional linear Kalman 
filter, there are well-known drawbacks (Julier and Uhlmann, 1997). The filter is prone to divergence if the 
underlying dynamical system is significantly nonlinear over the discrete time step, and the error 
covariance is often estimated poorly, which degrades performance and can cause divergence. Some 
researchers, including Julier and Uhlmann (1997), have developed improved algorithms similar to the 
EKF that are not limited by these issues. 


Estimating Bearing Properties in a Rotor-Bearing System 

The EKF will be applied to estimate the rotordynamic properties of two bearings in the rotor-bearing 
system depicted in figure 2. The bearings are assumed to have identical properties, which reduces the 
number of parameters to be estimated from sixteen to eight. As mentioned earlier, the geometry and 
operating parameters of the rotor-bearing system are chosen to be representative of NASA’s rotordynamic 
simulator test rig and are listed in table 1 , and the assumed stiffness and damping properties of the 
bearings listed in table 2 are thought to be realistic values for the chosen operating conditions. Shaft 
motion can be excited by the applied imbalance at the shaft ends or by impact or by both. In any case, the 
applied imbalance and the impact force are assumed to be known. The EKF algorithm will be tested by 
numerical simulation in which the “experimental” data for the shaft motion is generated by numerically 
solving the equations of motion, eq. (2), and then adding uncorrelated white noise with zero mean and 1.0 
pm standard deviation to simulate measurement noise. 


TABLE 1.— PARAMETERS OF THE ROTOR-BEARING SYSTEM 


h, h 

0.1245 m 

n 

20000 rp m 

d 1 , d 2 

0.1661 m 

O, r 2 

18.42 mm 

M 

3.089 kg 

m .\ , m 2 

0.5 g 

Ip 

2.629* 10 -3 kg-m 2 

<i>i 

0° 

It 

2.483* 10” 2 kg-m 2 

<t>2 

90° 
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TABLE 2.— STIFFNESS AND DAMPING PROPERTIES 


ASSUMED FOR THE BEARINGS 


k X x 

500 kN/m 

Cxx 

1,100 N-s/m 

k xv 

400 kN/m 

<Av 

50 N-s /m 

k\ X 

-350 kN/m 

c vx 

M0 N-s /m 

kyy 

600 kN/m 

Cyy 

1,200 N-s/m 


The performance of the EKF filter is largely affected by the values of user-defined tuning parameters 
and the initial estimates of the state and error covariance matrix. The Q and R matrices model the process 
noise and measurement noise, respectively. More process noise indicates to the filter that a greater 
difference exists between the filter model and the real process. Likewise, larger values in R indicate the 
presence of more error in the measurements. In general, these quantities can vary with time, but in this 
analysis, they are assumed to be constant. Their values are set by the user, and they are generally used as 
tuning parameters, especially for the EKF. The filter performance is also a function of the initial estimates 
of the state and error covariance matrix. Larger initial values of the error covariance matrix represent a 
greater uncertainty in the initial state estimates. 

In this application, the values of the process and noise matrices were tuned by a trial and error 
procedure. The EKF filter was generally insensitive to the values of the process noise for the 
rotordynamic coefficients, Q r , so it was set to zero, i.e., Q r = 0. Also, the filter worked best for zero 
values of the process noise for the state variables, i.e., Q v = Q z = 0, which is expected since the 
experimental data was generated with a model that is identical to the filter model. The trial and error 
procedure indicated that values for the diagonal elements of R between 0.2 and 2.0 pm produced 
acceptable results for the amount of measurement noise, so the diagonal elements of the measurement 
noise covariance matrix were set to R = diag [1.0 1.0 1.0 1.0] pm. 

The filter was much more sensitive to the initial estimates of the states and error covariance matrix 
than to Q and R. Initial estimates for the states were chosen to be 
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and the initial error covariance matrix was set to be a diagonal matrix with elements of 
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Results 

In each of the simulations that follow, the motion is simulated for 0.125 s, or approximately 41.7 
shaft revolutions, with a time step of 50 ps. The shaft originates from its equilibrium position with zero 
initial velocity. The first 30 ms of data made available for processing through the EKF is removed to 
eliminate the transients due to the initial conditions. 

For the first numerical simulation, the shaft motion was excited by both applied imbalance and 
impact. The imbalance is listed in Table 1, and the impact consisted of a constant force /,' = -515 N acting 
at the shaft mass center ( z p = 0), initiating at 40 ms and enduring for 50 ps. Figure 4 plots the percent 
error between the estimated coefficients and the actual values as a function of the time step for one 
representative simulation. Within 30 ms, the estimates for all but the cross-coupled damping reached 
steady state, with errors less than 0.5 percent for the direct damping and less than 5.5 percent for the 
stiffness terms. The errors for the cross-coupled damping varied much more significantly. This result is 
common among the parameter identification methods (De Santiago and San Andres, 2007a and 2007b). In 
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Time (ms) 

Figure 4. — Percent error in the estimated rotordynamic coefficients for a 
representative simulation. The shaft was excited by imbalance and impact. 



turbine end. Data shown includes the 
actual noisy “experimental” orbit and the 
orbit simulated with the rotordynamic 
properties estimated from the simulation 
with both imbalance and impact excitation 


general, the EKF algorithm improves with better initial guesses, so it is possible to use the EKF 
iteratively, especially if the estimates do not reach steady state in one filtering pass. In such a case, the 
final estimates from one pass through the filter would be used as the initial guesses for the next pass. This 
process is repeated until eventually the change in the final estimates is below some tolerance. In the 
present case, iteration did not improve the estimates significantly. 

Figure 5 shows a comparison between the orbit from a simulation using the estimated bearing 
properties and the actual noisy orbit data. One may be tempted to use this comparison as a check to verify 
the accuracy of the parameter estimates since the luxury of comparing the estimated bearing properties 
with the actual bearing properties will not generally be available. The logic is that a favorable orbit 
comparison indicates success of the identification scheme. In this case, the orbits match very closely, 
apparently indicating that the parameter identification was successful. The problem, however, is that the 
quality of the orbit prediction is not closely tied to the accuracy of each system parameter individually. 
Data presented later will show rather close orbits predicted using bearing parameters that are significantly 
in error compared to the actual values. 

Table 3 lists the results of the EKF applied to twenty numerical simulations where the shaft was 
excited by imbalance and impact. Considering the amount of noise in the data, the direct damping 
estimates were exceptional, the estimates of the four stiffness terms were reasonable, but the cross- 
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coupled damping estimates showed significant error. The deviation in the estimates (represented by the 
percent standard deviation relative to the actual parameter values) also followed the same trends. The 
cross-coupled damping estimates showed the largest relative variance, which is comparable to other 
estimation methods (Tiwari et al., 2004, De Santiago and San Andres, 2007b). 


TABLE 3.— AVERAGE ESTIMATED ROTORDYNAMIC COEFFICIENTS FROM 20 SIMULATIONS 


Parameter 

Unit 

Exact 

value 

EKF 

estimate 

Standard 

deviation 

Percent 

error 

Percent standard 
deviation 

Cxx 

N-s/m 

1100 

1095.0 

9.2 

-0.455 

0.840 

<Av 

N-s/m 

50 

39.6 

9.3 

-20.8 

18.6 


N-s/m 

^to 

-73.4 

11.2 

83.5 

28.0 

c w 

N-s/m 

1200 

1208.3 

13.5 

0.692 

1.13 

kxx 

N/m 

500,000 

520,094 

21,644 

4.02 

4.33 

kxv 

N/m 

400,000 

388,402 

20,732 

-2.90 

5.18 

kyx 

N/m 

-350,000 

-311,207 

21,843 

-11.1 

7.28 

kyy 

N/m 

600,000 

533,420 

21,110 

-11.1 

3.52 


The results shown here indicate that the EKF is capable of estimating the rotordynamic coefficients 
when the shaft motion results from imbalance and impact excitation, but conditions do not always exist 
where sufficient access to the shaft is available to deliver an impact to the shaft. In such cases, it would be 
advantageous to estimate the bearing properties with data from imbalance excitation only, for it is more 
common to have knowledge of the shaft residual imbalance even if it is not possible to apply a known 
imbalance. In this effort, the EKF filter was applied to identify the bearing properties when the system 
was simulated with imbalance excitation only, and the errors in the EKF estimations are shown in figure 
6. Even though the errors in the direct stiffness estimates were reasonably small (less than 8 percent), the 
errors in the cross-coupled stiffness estimates were considerably larger, ranging from 24 to 50 percent, 
and the errors in the cross-coupled damping estimates were much worse (>100 percent). Even after 
significant effort to tune the user-definable parameters, the estimates from the EKF filter could not be 
materially improved. The conclusion is that the EKF filter cannot estimate the bearing properties without 
the impact excitation present. Some parameter identification methods (Lee and Hong, 1989; Tieu and 
Qiu, 1994; Tiwari et al., 2002; De Santiago and San Andres, 2007b) are capable of estimating bearing 
properties using imbalance excitation only, though the results are not as accurate as when impact 
excitation data is used. It is possible that this particular rotor-bearing problem exacerbates the inherent 
drawbacks with the EKF and that other Kalman filtering methods that overcome these problems, such as 
the Unscented Kalman Filter developed by Julier and Uhlmann (1997), may be more successful at 
identifying bearing parameters from imbalance response data. 

It is interesting to compare the orbits computed from the (significantly incorrect!) estimated 
parameters and the true, noisy experimental data. One may expect that the two orbits would be noticeably 
different since the estimated parameters are significantly in error, but figure 7 shows that the two orbits 
are, instead, very close, as close as the orbits with the much more accurate bearing parameter estimates 
(see fig. 5). Even though the EKF cannot accurately distinguish the bearing properties individually in this 
case, the filter still estimates the bearing parameters collectively in such a way that the net effect of the 
estimated bearing properties on the shaft orbit is similar to the true effect of the actual bearing properties. 
Consequently, it is not advisable to use orbit comparison as a method to assess the quality of the 
parameter identification method. 


NASA/TM— 2008-215298 


10 




300 


200 
•„ 100 
I 01 


-100 
-200 1 


-300 



^ 555 §=:: 


M 


30 


Figure 6- 


“ I - 

45 


—r~ 

60 


"T" 

75 


— CV™ Cxx 

--O-- k-xx 

— o— Cxy 

--□--kxy 

Cyx 

--A— kyx 

— O-Cyy 

"O-Kr 


90 


r 

105 


120 


Time (ms) 

-Percent error in the estimated rotordynamic coefficients. 
The shaft was excited by imbalance only. 



turbine end. Data shown includes the 
actual noisy “experimental” orbit and the 
orbit simulated with the rotordynamic 
properties estimated from the simulation 
with only imbalance excitation. 


Conclusions 

This paper describes an Extended Kalman Filter (EKF) designed to identify stiffness and damping 
properties of two identical journal bearings from measurements of the shaft motion in a rotor-bearing 
system. The EKF successfully identified the rotordynamic bearing properties in numerical experiments 
when the shaft motion was excited by impact and imbalance, but it failed when the shaft motion was 
excited solely by imbalance. The filter performance was influenced marginally by user-defined tuning 
parameters, Q and R, and more significantly by the initial values of the bearing parameters and error 
covariance matrix. Accuracy of the parameter estimates was similar to other published bearing parameter 
identification methods. Estimates for the direct damping and direct and cross-coupled stiffness values 
were considerably better than estimates for the cross-coupled damping values. Likewise, the relative 
variance in the cross-coupled damping estimates over several simulations was considerably larger than 
the relative variance in the direct damping and direct and cross-coupled stiffness estimates. 

The paper also discusses the inadequacy of orbit comparison as a means to assess the quality of 
bearing parameter estimates from the EKF. When the EKF converged to a set of bearing parameter 
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estimates, the shaft orbits determined analytically using those parameter estimates were close to the actual 
shaft orbits, regardless of the quality of the individual bearing parameter estimates. 
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Appendix — Symbols 
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r 0 
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s k 


linearization of f about the current state estimate 

system damping matrix 

total applied forces on the shaft 

system gyroscopic matrix 

output matrix 

identity matrix 

transverse, polar shaft mass moment of inertia 
Kalman gain matrix 
length between the bearing centers 
mass matrix 

continuous-time estimate of the error covariance matrix 
estimate of the error covariance matrix at time step k 
initial estimate of the error covariance matrix 
prediction of the error covariance matrix at time step k 

covariance of the noise vectors corresponding to the shaft velocity, shaft displacement, 
and bearing properties 

covariance of the global system process noise vector 

covariance of the measurement noise vector 

system stiffness matrix 

linearized bearing damping coefficients 

distance from the shaft center of mass to balance plane 1 , 2 

nonlinear function representing the total system dynamics 

x, y direction bearing force 

x, y direction impact force 

nonlinear function representing the shaft dynamics 
linearized bearing stiffness coefficients 

distance from the shaft center of mass to the centers of the bearings 
shaft mass 

total imbalance mass (residual and applied) at balance plane 1 , 2 

distance from the shaft center of mass to the displacement probes 

eccentricity of the imbalances at balance plane 1, 2 

state vector of rotordynamic coefficients 

initial estimates of the rotordynamic coefficients 

global state vector 

global state vector at time step k 
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s continuous-time estimate of the global state vector 

s /( estimate of the global state vector at time step k 

s 0 initial estimate of the global state vector 

s* prediction of the global state vector at time step k 

s v , s z state vector of velocity and displacement variables 

t time 

t k time at time step k 

y k output measurement vector at time step k 

v k measurement noise vector at time step k 

w v , w z , w r noise vectors corresponding to the shaft velocity, shaft displacement, and bearing 
properties 

w global system process noise vector 

X\, y\ x, y shaft displacement at the bearing center near the turbine disk 

x 2 , >’2 x, y shaft displacement at the bearing center near the thrust bearing plate 

Jtio, >’ i o equilibrium x, y coordinate of the shaft at the probe near the turbine disk 

X20, Via equilibrium x, y coordinate of the shaft at the probe near the thrust bearing plate 

x p 1, y p \ x, y shaft displacement at the probe near the turbine disk 

Xp 2 , y p 2 x, y shaft displacement at the probe near the thrust bearing plate 

z,z,z state vector of shaft displacement, velocity, and acceleration variables 

z 0 initial value of the shaft displacement state vector 

z 0 initial value of the shaft velocity state vector 

z p distance from the shaft center of mass to the location of impact excitation 

4*1, <(>2 imbalance angle measured from the reference axis at balance plane 1 , 2 

Q shaft speed 
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